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ABSTRACT 

In this work, we present a detailed study of the rotational properties of mag- 
netized and self-gravitating dense molecular cloud cores formed in a set of two 
very high resolution three-dimensional molecular cloud simulations with decay- 
ing turbulence. The simulations have been performed using the adaptative mesh 
refinement code RAMSES with an effective resolution of 4096^ grid cells. One 
simulation represents a mildly magnetically-supercritical cloud and the other a 
strongly magnetically-supercritical cloud. We identify dense cores at a number 
of selected epochs in the simulations at two density thresholds which roughly 
mimick the excitation densities of the NH3 (J — K) = {1,1) transition and the 
N2H+ (1-0) emission line. A noticeable global difference between the two simu- 
lations is the core formation efficiency (CFE) of the high density cores. In the 
strongly supercritical simulations the CFE is 33 percent per unit free-fall time of 
the cloud (tff^ci), whereas in the mildly supercritical simulations this value goes 
down to ~ 6 percent per unit A comparison of the intrinsic specific angular 

momentum (jsz?) distributions of the cores with the specific angular momentum 
derived using synthetic two-dimensional velocity maps of the cores (j2Z)), shows 
that the synthetic observations tend to overestimate the true value of the specific 
angular momentum by a factor of ~ 8 — 10. We find that the distribution of 
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the ratio j3D/j2D of the cores peaks at around ~ 0.1. The origin of this discrep- 
ancy hes in the fact that contrary to the intrinsic determination of j which sums 
up the individual gas parcels contributions to the angular momentum, the de- 
termination of the specific angular momentum using the standard observational 
procedure which is based on a measurement on the global velocity gradient under 
the hypothesis of uniform rotation smoothes out the complex fiuctuations present 
in the three-dimensional velocity field. Our results may well provide a natural 
explanation for the discrepancy by a factor ~ 10 observed between the intrinsic 
three-dimensional distributions of the specific angular momentum and the corre- 
sponding distributions derived in real observations. We suggest that previous and 
future measurements of the specific angular momentum of dense cores which arc 
based on the measurement of the observed global velocity gradients may need to 
be reduced by a factor of ~ 10 in order to derive a more accurate estimate of the 
true specific angular momentum in the cores. We also show that the exponent of 
the size-specific angular momentum relation are smaller (~ 1.4) in the synthetic 
observations than their values derived in the three-dimensional space 1.8). 

Subject headings: ISM : clouds - ISM : globules - ISM : kinematics and dynamics 
- ISM : magnetic fields - turbulence - MHD 



1. INTRODUCTION 

An important issue for current ideas of star formation is whether most stars of a given 
mass are born single, in a binary system or in multiple systems (e.g., Bodenheimer 1995; 
Bodenheimer et al. 2000; Larson 2010). Among other several physical processes that affect 
the fragmentation of cores prone to star formation such as their geometry (e.g. Boss 2009), 
their metallicity (e.g., Hocuk & Spaans 2010), and degree of magnetization (e.g.. Boss 1999; 
Price & Bate 2007; Hennebelle & Teyssier 2008) the rotational properties of the cores will 
also strongly affect their ability to fragment (e.g.. Nelson 1998; Sigalotti & Klapp 2001; 
Matsumoto & Hanawa 2003; Hennebelle et al. 2004; Machida et al. 2009; Walch et al. 2009). 
This may have important consequences on the fraction of binary and multiple systems. The 
determination of the angular momentum of the cores and the fraction of their energy stored 
in rotational motions is complicated by the fact that the velocity fields in the cores usually 
exhibit complex supersonic motions when mapped with tracers such as CO, CS, and C^^O 
with a transition to subsonic, coherent motions when tracers with a higher excitation density 
are used such as NH3 and N2H+ lines (e.g. Barranco & Goodman 1998; Goodman et al. 
1998; Caselh et al. 2002a). Usually, the assumption is often made in the observations that 



-3- 



cores have a uniform rotation and follow a rigid-body rotation law. Their angular velocity 
Q is deduced from measuring global velocity gradients in velocity maps built using velocity 
measurements along the line of sight. Such measurements for a large number of molecular 
clouds, clumps or cores were performed by Myers & Benson (1983), Goldsmith & Arquilla 
(1986), Kane & Clemens (1997), Pound & Goodman (1997), and Rosolowsky (2007) using 
CO lines, and by Goodman et al. (1993), Barranco & Godman (1998) using observations 
in the (J - K) = {1-1) transition of NH3, and CaseUi et al. (2002a), Pirogov et al. (2003), 
Olmi et al. (2005) and Chen et al. (2007) using N2H+ (1-0) hne observations. Several other 
authors have also obtained rotation measurements for individual cores or small numbers of 
cores (e.g., Harris et al. 1983; Menten et al. 1984; Armstrong et al. 1985; Wadiak et al. 
1985; Zheng et al. 1985; Ho & Haschick 1986; Jackson et al. 1988; Ho ct al. 1994; Ohashi 
et al. 1997; Bcllochc et al. 2002; Caselli et al. 2002b; Tafalla et al. 2004; Shinnaga et al. 
2004; Redman et al. 2004, Schnee et al. 2007; Chen et al. 2008, Chen et al. 2009; Csengeri 
et al. 2010). 

It is important to investigate whether the distributions of the rotational properties of 
cores show significant variations as a function of the adopted density tracer /threshold and 
as a function of the environment (e.g., for example for different magnetization levels of the 
parent cloud), and most importantly whether the distributions that are determined from ob- 
servations are a faithful reproduction of the distributions in the intrinsic three-dimensional 
space, would the entire dynamical and structural properties of the cores be accessible. Fur- 
thermore, knowing the statistical distributions of the rotational properties of dense cores in 
molecular clouds is crucial in order to asses the statistical relevance of rotational parameters 
assigned to individual rotating core collapse simulations. In previous works, Burkert & Bo- 
denhcimer (2000) superimposed a random velocity field of power spectrum P{k) oc A;" (with 
n between —3 and —4) onto the density field of turbulent molecular cloud cores. They found 
that the projected velocity maps display velocity gradients that can be interpreted as rotation 
and measured the average specific angular momentum of the cores to be 7 x 10^'^(i?c/0.1)^'^ 
cm^ s"^, where Rc is the size of the core, and the average value of the rotational parameter 
(i.e., ratio of rotational energy to gravitational energy) of /3rot ~ 0.03. Gammie et al. (2003) 
examined the distribution of specific angular momentum for cores formed in a set of magne- 
tized, self-gravitating molecular clouds with decaying turbulence at the resolution of 256^. 
They showed the distribution of specific angular momentum at an intermediate epoch in one 
of their simulation (the one with beta plasma value of f3p = 0.1 after 0.09 sound crossing 
time). The specific angular momentum they found follows a distribution with values ranging 
between 10^^ cm^ s~^ to 7 x 10^"^ cm^ s^^ with the peak of the distribution at j = 4 x 10^^ 
cm^ s~^. Jappsen & Klessen (2004) found that there might be a dependence of the average 
specific angular momentum of the cores as a function of the Mach number of the cloud. Li 
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et al. (2004) and Offner et al. (2008) also evaluated the specific angular momentum j and 
rotational parameter Prot for cores formed in their molecular cloud simulations, both with 
resolutions of 512^ grid cells, and found values of the order of the observed ones. 

We revisit the issue of rotation in molecular clouds cores using very high resolution 
simulations of magnetized, turbulent, and self-gravitating isothermal molecular clouds (MCs) 
performed with the Adaptative Mesh Refinment code RAMSES. Our primary aim is to 
quantify any systematic difference between the intrinsic angular momentum of the cores and 
the values of the specific angular momentum derived from synthetic observations built using 
an approach that mimicks the procedure commonly used in the observations. In § [2|, we 
describe the molecular cloud simulations and the clump finding algorithm used to identify 
dense cores in the simulations at different density thresholds. In §[3] we describe some of the 
basic properties of the simulations. The rotational properties of dense cores are presented 
and discussed in § H] and in § [5] we summarize our results and conclude. 



2. SIMULATIONS AND CLUMP-FINDING 

In this study, we analyze two three-dimensional, isothermal, magnetized, and self- 
gravitating molecular cloud simulations with decaying turbulence. The simulations were 
performed using the adaptative mesh refinement (AMR) code RAMSES (Teyssier 2002) on 
a cubic grid with periodic boundary conditions in all directions. RAMSES is a second order 
Godunov scheme and uses the constraint transport method to ensure divB = (Fromang 
et al. 2006). The initial setup of the simulations is similar to those presented by Dib et 
al. (2007a,2007b, 2008a, 2008b) on uniform grids with 256^ and 512'^ resolutions. However, 
by using the AMR technique, the linear resolution is enhanced by almost an order of mag- 
nitude which is extremely important in order to study the properties of dense cores which 
occupy but a small fraction of the simulation volume. The initial coarse grid is uniform 
and consists of 512^ cells. The density field on the initial grid is homogeneous. Refinement 
of the grid was performed such as to keep a local Jeans mass resolved with 10 grid cells. 
The grids were refined for up to three additional levels leading to an effective resolution of 
4096^ grid cellfl In contrast to our previous simulations in which turbulence was contin- 
uously forced, in the current simulations, turbulent motions that are injected at the start 
of the simulations are left to decay freely as time evolves. The initial local components of 



-'^However, at the last time-steps included in our analysis, the number of grid cells in the most refined 
level remains quite modest and is of the order of a few thousand cells (e.g., 5311 cells at the last considered 
time-step in simulation B2). 
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the turbulent velocity field are modeled using a fractal brownian motion power spectrum 
with an exponent of -11/3, corresponding to a Kolmogorov spectrum. The initial velocity 
perturbations are generated on all scales of a 256^ cells grid and interpolations are performed 
to acquire the velocity components on the initial coarse 512^ grid. The amplitudes of the 
velocity components are chosen such as to represent an initial rms sonic Mach numbei@ of 
Ai = 10. The Jeans number of the box is Jbox = 4 (i.e., number of Jeans masses in the box 

is Mhox/Mjeans,box = Jbox = 64). 

In physical units, the simulation box has a size of L^i = 4 pc (thus the spatial resolution 
is 4 pc/4096 ~ 10~^ pc), an average number density of Uaver = 500 cm~^, and a temperature 
of T = 11.43 K which leads to a sound speed of 0.2 km s~^ with a mean molecular weight of 
Myj = 2.36 in units of the proton mass ttih- The free-fall time of the simulation volume is 
defined as being tjf^d = 1/(G' pY^"^ ~ 2.76 x 10^ yrs, where p is the physical uniform density 
in the bojjf]. The simulations differ by the strength of the initial magnetic field in the parent 
cloud, with one simulation being slightly magnetically-supercritical (run Bl), and the other 
one strongly supercritical (run B2). Initially the magnetic field has a component only in one 
direction. The initial magnetic field strength, the mass-to-magnetic flux ratio (in units of 
the critical value for collapse fXcr ~ (4 vr^ G)~^/^, Nakano & Nakamura (1978 )@ and the beta 
plasma values are = 14.5,4.6 /iG, = 2.8,8.8, and (3p = 0.1, 1, for simulations Bl, and 
B2 respectively. The Poisson equation is solved, right from the beginning of both simulations, 
using the conjugate gradient method implemented in RAMSES. In this work, we restrict our 
analysis to magnetized cloud models since, aside from the now commonly accepted fact that 
clouds and cores are indeed magnetized (e.g., Crutcher et al. 2004), we have shown in some 
of our previous work that cores formed in magnetized cloud provide the best match to the 
available observational constraints such as the exponent of virial parameter-mass relation of 
the cores (Dib et al. 2007a) and their lifetimes (Galvan-Madrid et al. 2007) (particularly 
for the nearly magnetically-critical cloud). 

In order to be able to handle the large data cubes when performing the clump finding. 



^The initial Alfven Mach number in simulations Bl and B2 is thus equal to M-a — A4 x j3p — 0.1 and 
10, respectively 

■^The value of i//.cz defined here is larger than the free-fall time of a sphere of uniform density p by the 
factor (32/3 tt^^^ 

^As in our previous work, we use the value derived by Nakano & Nakamura (1978) for ficr, which is 
derived from a linear stability analysis and which is potentially more adapted to our initially uniform density 
configuration threaded by a unidirectional magnetic field. The value of ficr = 1/(63 G)^/^ derived by 
(Mouschovias & Spitzer 1976) from non-linear calculations differs from the first quantity by a factor of 25 
per cent. 
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we have binned the simulations such as to generate 1024^ uniform data cubes for each of 
the time-steps of the simulations we have analyzed. Cores are identified on the uniform 
grids using a clump-finding algorithm based on a density threshold and a friend-of-friend 
criterion similar to the ones used in Dib et al. (2007,2008a,2008b), Dib & Kim (2007), Audit 
& Hennebelle (2005), and Hennebelle & Audit (2007). We have performed a clump finding 
directly in the three-dimensional space using a variety of density thresholds ranging from 
10 Haver to 1280 Haver with increments by a factor of 2. We make the simple assumption 
that a threshold of 20 Haver = 10^ cm~^ is of the order of the excitation density for the NH3 
(J — K) = {1,1) transition, whereas a threshold of 160 Haver = 8 x 10^ cm~^ is of the order 
of the excitation density for the (1-0) N2H"'" emission line. Below these values, it is assumed 
that the molecules are under-excited and that there is no emission in the considered lines. 
The clump-finding algorithm takes into account the periodicity of the simulation box and 
thus, of cores that extend across the box boundaries. 



3. GENERAL PROPERTIES OF THE SIMULATIONS 

In this section, we briefly summarize some of the basic properties of simulations Bl 
and B2 in terms of their respective populations of cores identified at the selected den- 
sity thresholds. As time evolves, turbulence decays in the box and the sonic Mach num- 
ber decreases from = 10 at t = to M = [5.28, 3.54, 2.97, 3.08, 3.28, 3.36] at t = 
[0.087, 0.226, 0.293, 0.356, 0.410, 0.424] tfj^^i is run Bl and to A^ = [3.61, 3.37, 2.89, 2.89, 2.59, 
2.73] at t = [0.257,0.280,0.351,0.373,0.420,0.486] in model B2. By t ~ 0.42 iff^^i, the rms 
Mach number has decayed faster in model B2 to a value of 2.6 and which is comparable to the 
rms Mach number in the Ophiucus cloud (e.g., Enoch et al. 2007 based on the COMPLETE 
survey. Ridge et al. 2006) as compared to a value of ~ 3 in run Bl which is comparable 
to the rms Mach number in the Perseus molcular cloud (Enoch et al. 2007, Pineda et al. 
2008)@. In both runs, the Mach numbers increases again at later stages because the rms 
velocity becomes increasingly dominated by the large velocity components in the collapsing 
cores. Tabs. [T]|l] display the rms Mach numbers, the numbers of C20 cores (cores identified 
at a density threshold of 20 Haver) and the C160 cores (cores identified at a density threshold 
of 160 Haver), the peak density in the box in units of the average density, the median mass of 
the cores, their median sizes (taken to be the cubic root of the cores volumes), their median 



■^The velocity dispersions quoted by Enoch et al. (2007) in the Perseus and Ophiucus molecular clouds 
are averaged values measured from all the pointings in each cloud. This averaged velocity dispersion, and the 
associated Mach number calculated assuming a gas temperature of 10 K can be associated with the average 
'width' of the cloud along the line of sight which is expected to be of the order of a few pc 
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specific angular momentum, and their median rotational parameter (see below for definition) 
at two selected, nearly similar, epochs. At early epochs (i.e., when t < 0.35 tjf^ci), the peak 
density in the simulation box, associated with the most massive cores is of the order of a few 
10^ Uaver ~ a few 10^ cm~^ which is characteristic of the densities of gravitationally bound 
cores (e.g., Lee et al. 2001; Kandori et al. 2005). At t ~ 0.42 the most massive cores in 
the simulations are in a stage of free-fall collapse and the peak density is of the order of a few 
10^ Uaver ~ ci fcw 10'' cm~'^ which is similar to the densities of cores showing signature of infall 
(e.g., J0rgensen et al. 2002). Fig. [T] shows the time evolution of the number of cores in the 
simulation box in models Bl and B2 at the thresholds of rithr = 10,20,40,80, 160,320,640, 
and 1280 Uaver, respectiveljl^. Fig. [2] displays the time evolution of the mass fraction of the 
box enclosed in the C20 cores and the C160 cores in the two simulations. Fig. [3] and Fig. H] 
display the column density, along the three directions of the box, of a few selected C160 
cores at t = 1.152 Myrs = 0.417 tff^^i in simulation Bl and at t = 1.159 Myrs = 0.420 tfj^^i 
in simulation B2, respectively. The physical size of the cores is indicated by a marking of the 
0.01 pc scale at the lower left corner of each map. The cores show a variety of morphologies 
going from more roundish structures to more filamentary ones. Some of the cores exhibit 
a single density peak whereas others are already fragmented into two or several sub-cores. 
Several authors (e.g.. Price & Bate 2007, Hennebelle & Teyssier 2008) have suggested that 
strong perturbations must be seeded in the clumps/cores in order to allow them to fragment 
in the presence of a magnetic field. Our simulations suggest that such strong perturbations 
which would allow for fragmentation, occur naturally in cores formed in a gravo-turbulent, 
magnetized, molecular cloud. 

There are a number of similarities between the two simulations and a number of no- 
ticeable differences. In both simulations, the number of cores identified at the lower density 
thresholds {rithr = 10 — 80 Uaver) declines as time advances, while the number of cores iden- 
tified at the higher density threshold {rithr = 160 — 1280 Uaver) increases. The increase in 
the high density cores is simply due to the effect of gravity which generates more centrally 
peaked objects as time advances. On the other hand, the decrease in the low density cores, 
in addition to being a signature that some of them are not gravitationally bound and may 
disperse as earlier suggested by Vazquez-Semadeni et al. (2005a) and Dib et al. (2007a), is 
also an indication of core coalescence which depletes a fraction of these low density cores (Dib 
et al. 2007b, Dib 2007, Dib et al. 2008c). The number of cores in the stronger magnetized 



^The data starts in figures 1, 2, and 3 starts at slightly different points in time as we have dumped the 
output at equal number of the simulations time-step which varies in time and from simulation to simulation. 
The data points were chosen such as to cover the epoch at which dense structures form and up to the point 
where a runaway collapse starts to set in in one of the cores. 
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cloud simulation (Bl) is observed to be initially smaller than in the strongly magnetically- 
critial cloud simulation B2 and the number of dense cores formed at the later stages in the 
evolution of the cloud is also observed to be smaller in the Bl cloud simulation than in the 
B2 model. This is also consistent with previous findings (Vazquez- Semadeni et al. 2005a 
and Dib et al. 2007a) which suggests that the number of cores formed in the cloud decreases 
when the cloud is more magnetically supported. 

As can be observed in Fig. [21 the mass fraction of the cores with respect to the mass 
enclosed in the box, e = Mcores/ M^ox is larger and increases at a faster pace in the strongly 
supercritical simulation B2 than in the mildly supercritical one, Bl. The definition of the 
core formation efficiency as being the mass enclosed in the cores per unit free-fall time of 
the cloud: 

CFE = (1) 

yields core formation efficiencies for the C20 and C160 cores, based on the variations of e 
between the first and last time-step of each simulation in Fig. [21 of CFE (C20) = 0.29 and 
CFE (C160) = 0.059, and CFE (C20)= 0.41 and CFE (C160) = 0.33 for simulations Bl 
and B2, respectively. The results for the strongly supercritical simulation (B2) are consistent 
with previous results (e.g., Li & Nakamura 2004; Vazquez-Semadeni et al. 2005b, Nakamura 
& Li 2005). Our results also show that the CFE is strongly reduced with an increasingly 
important magnetic field in the cloud, particularly for the case of the higher density cores. 
It is useful to emphasize any comparisons between isothermal models such as ours and the 
observations remain valid as long as we are dealing with pre-stellar cores which are themselves 
embedded in a cloud that does not harbor strong sources of radiation (stars/protostars). The 
presence of radiation will unavoidably alter the structure of the cloud and reduce the core 
formation efficiency, and subsequently the star formation efficiency as recently pointed out 
by the simulations of Bate (2009) and Offner et al. (2009). 



4. ROTATIONAL PROPERTIES OF MAGNETIZED CORES 

After having defined the cells that belong to the same core, we calculate the components 
of the angular momentum of the core, Jc around the center of mass using the standard 
expression: 

Jc = ^ X mj Vj, (2) 
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where r is the distance from each cell to the center of mass, rui the mass of the cell, and 
Vj is the relative velocity of the cell with respect to the core's center of mass. The specific 
angular momentum is then calculated as being: 



j^n = ^, (3) 



where Mc is the mass of the core. The rotational parameter ^^ot, calculated in the 
three-dimensional space, is usually defined as being: 



%ot = (4) 

I ^grav \ 



where Egrav and Erot are the core's gravitational energy and kinetic energy stored in 
rotational motions, respectively. The latter is given by its basic definition: 



Erot^ I IprWdV, (5) 

'Vc ^ 



where Vc is the volume of the core, and io is the local angular velocity given by: 



The gravitational energy Eg^av of the cores is observationally approximated by Egrav = 
—qGMl/Rci where Mc is the mass of the core and is a positive number which accounts 
for the core morphology and inner mass distribution. However, as ^^ot is intended to rep- 
resent a volume integrated estimate of the balance between centrifugal forces and effective 
gravitational forces taking into account that dense cores in molecular clouds are not isolated 
objects, and in order not to make any assumption on the values of q in the 3D estimate, it is 
more adequate, from a theoretical point of view, to replace Eg^-av by the quantity W which 
is slightly different than Eg^av (e.g., Dib et al. 2007a) and which is given by: 



where (j) is the gravitational potential. Note that W is not exactly equal to the volume 
gravitational energy. -Egrac- because the true gravitational potential is a result of the dis- 
tribution of matter inside the core and outside of it (e.g., the parent cloud). However, as 
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already shown in Dib et al. (2007a), the essential part of the gravitational acceleration in 
the core is due to the mass contained inside the core and in general W ~ Eg^av Thus, our 
definition of Prot is: 

Prot = Erot/\Wgrav\- (8) 

It is important to stress that in a time dependent model of a turbulent cloud with gravity, 
the statistical properties of the cores formed in the cloud may vary with time (e.g., Jappsen 
& Klessen 2004; Dib et al. 2007a). In addition to angular momentum loss that may be 
due to the cores secular evolution (i.e., magnetic braking), the angular momentum of a core 
might be modified if the core enters a soft gravitational encounter with another/other core(s) 
and/or a physical merger (e.g., Larson 2010). Following a collision, the remnant merger may 
inherit a reduced or enhanced specific angular momentum. In our simulations. Fig. [1] shows 
that the number of cores detected at the lowest thresholds declines as a function of time 
while the number of cores identified at the highest thresholds increases. As some of the cores 
evolve into a stage of gravitational contraction, their average density will increase and they 
will be detected at higher density thresholds. However, as most of the cores identified at the 
thresholds of 10 — 80 Uaver are gravitationally unbound, and have generally long dispersion 
timescales (i.e., Vazquez-Semadeni et al. 2005), it is likely that the rapid decline in their 
numbers is due to merger events. Thus, in a dynamically evolving cloud, it is important to 
keep in mind that the overall angular momentum of the dense cores that form in the cloud 
will be re-distributed among the cores and the inter-core medium. 

4.1. Distributions of intrinsic specific angular momentum and rotational 

parameter 

Figs. [5] and |6] display the time evolution of the normalized distributions of the intrinsic 
specific angular momentum j^D and of the rotational parameter fSrot for the C20 and C160 
cores identified in the two simulations. Albeit our models were not necessarily tailored to 
mimick any particular star forming region, it is temptative to over-plot to the numerical 
results the observationally derived NH3 and N2H+ distributions of j and /3rot (dashed-line 
histograms). The NH3 cores are those observed by Goodman et al. (1993) and Barranco & 
Goodman (1998), whereas the N2H"'" cores are those observed by Caselli et al. (2OO2a)0. The 
median value of j^c for the C20 and C160 cores fiuctuates in time but remains close to jsD.med 



^It should be noted that the NH3 and N2H+ cores of Goodman et al. (1993) and Caselli et al. (2002a) 
and which are drawn from the sample of cores built by Benson & Myers (1989) contain a mix of cores from 
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(C20)=2 - 3 X 10^0 cm2 s"^ and j3D,med (C160)=8 x 10^^ - 1.5 x 10^° cm^ s'l in simulation 
Bl and j3D,med (C20)=2 - 3 x 10^° cm^ and j3D,med (C160) = 1.5 - 2 x 10^° cm^ s'^ in 
simulation B2. Similarly the median value of Prot for the two populations of cores fluctuates 
in time in the range of Prot,med (C20) = (1.8-3) x lO'^ and Prot,med (C160) = (0.48-0.8) x lO'^ 
in simulation Bl and (3rot,med (C20) = (1.4-2.3) x lO'^ and f3rot,med (C160) = (0.34-0.68) x 10"^ 
in simulation B2. Fig. [7] displays the time evolution of the median value of the C20 and C160 
cores Jsd's in simulations Bl and B2, and Fig. |8] displays the time evolution of the median 
values of Prot for the two populations of cores in the two simulations. The typical median 
absolute deviations of jso, shown as error bars in Fig. [71 are of the order of 40-50 percent 
and reflect the large width of the distributions in Figs. |5l The median absolute deviations of 
Prot, shown as error bars in Fig. |8l are of the order of 30-40 percent. The median values of 
the specific angular momentum and of the rotational parameter for the C20 and C160 cores 
in the two simulations Bl and B2 are compared to those of the observed NH3 cores and 
N2H^ cores in Fig. [7] and Fig. |8l respectively. Albeit there is a certain overlap in the median 
values of jso between simulations Bl and B2, the fact that the cores (particularly the C160 
cores) in the more magnetized model have smaller values of jsD is consistent with the fact 
that magnetic braking leading to a loss of angular momentum is playing a more important 
role in the more strongly magnetized cloud model (Bl) than in the less strongly magnetized 
cloud model (B2). We will address the issue of the origin of angular momentum loss by 
evaluating the effects of magnetic braking and the effects of gas accretion and coalescence of 
cores in a separate, dedicated, paper. 

A striking point in the comparison between the numerical specific angular momentum 
distributions of the C20 and C160 cores and the observed NH3 and N2H+ distributions 
(i.e.. Fig. |5] and Fig. [7]), respectively, is that the observations appear to overestimate the 
amount of the specific angular momentum by a factor of < 10. A similar trend is observed 
in the Prot distributions in Fig. IH] and FigJH] albeit the shift in this case appears to be of 
the order of a factor ~ 2. Since, the observational determination of Prot relies on the 
determination of the specific angular momentum of the cores but also on the determination of 
the gravitational energy, many uncertainties can affect the determinations of Prot such as the 
assumed idealized morphologies of the cores and the idealized distribution of matter within 
them. In the following, we will therefore focus our attention on the discrepancy between 
the observed and numerically measured specific angular momentum distributions. Offner et 
al. (2008) calculated the values of the specific angular momentum for some of the cores in 
their simulations (at the resolution of 512^) both in 3D and from projected velocity maps 



different star forming regions including Ophiucus and Taurus that are located at different distances and 
which possess different dynamical characteristic (e.g., different Mach numbers). 
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of the cores. They found that the 2D cores have specific angular momenta values that are 
systematically larger by a factor of ~ 10 than the values of j for the 3D cores. Several effects 
may induce differences between a numerically calculated distribution of specific angular 
momentum and an observed one. In the present simulations, we used an initial rms Ma 
number of 10 and which decays to ~ 3 by t ~ 0.42 tff^d whereas the observed cores 
by Goodman et al. (1993) and Caselli et al (2002a) are sampled from different molecular 
cloud complexes in which the dynamical conditions (e.g., the Mach number) are different. 
The most simple situation is that the parameters governing the dynamical evolution of the 
molecular cloud in the models do not exactly match the dynamical conditions in the regions 
where the observed cores are formed and thus cores in the observations may have intrinsically 
more specific angular momentum than the simulated ones. The second effect may be related 
to the fact that the simulated cores are identified in intrinsic space whereas the observed 
cores are identified in position-position-velocity cubes. Although some differences in the 
cores populations may arise when using the two methods, most of the most massive cores 
will be detected in the PPV method with their masses modified by a correction factor of order 
unity (e.g.. Smith et al. 2008). It is therefore unlikely that a factor of ~ 10 difference can 
originate from variations from the clump finding algorithm. A third possibility is that the 
totally different ways of calculating the specific angular momentum in the observations and 
in the simulations may stand behind this discrepancy. In order to test the latter hypothesis, 
it is necessary to generate synthetic velocity maps of the cores using individual projections, 
such as to eliminate any potential effect of blending of the cores along the line of sight, 
and measure the specific angular momentum from the projected velocity maps following the 
standard observational procedure. Unlike previous simulations, our models which have a high 
spatial resolution, and a large number of cores allow us to compare the entire distributions 
of the intrinsically measured specific angular momentum distributions to the ones derived 
from synthetic observations of the cores. 

We generate synthetic velocity maps of the C160 cores in both simulations at a nearly 
similar epoch (i.e., at t = 1.152 Myrs = 0.417 t//,d for simulation Bl and t = 1.159 Myrs 
= 0.420 for simulation B2) along the three main directions of the box. As stated 
above, the cores are projected individually such as to eliminate any effect of blending along 
the line of sight. The velocity in each pixel of the cores is the mean velocity in the line 
of sight. The velocity maps corresponding to the cores shown in Fig. [3] and Fig. H] are 
displayed in Fig. [9] and Fig. [TOl respectively. The velocity maps exhibit a variety of features 
ranging from well ordered motions that can be assimilated to rotation, both for roundish 
and filamentary cores, to more complex dynamical configurations in which there is no clearly 
verifiable velocity gradient. In general, the velocity maps of the cores, especially when the 
same cores are seen along the three different projections, do not support the idea of simple 
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rigid-body rotation. It is interesting to note that recent high spatial and spectral resolution 
observations using the Plateau de Bures interferometer by Csengeri et al. (2010) of five 
massive dense cores in Cygnus-X show a variety and level of complexity in their dynamical 
features similar to the ones observed in our synthetic velocity maps. 

In order to derive the velocity gradient from the synthetic velocity maps, we employ the 
program VFIT, developed by Goodman et al. (1993), which calculates the global velocity 
gradient in the map, Q and its orientation. The uncertainty in the centroid velocity needed 
as an input by VFIT is taken to be the velocity dispersion of the velocity distribution along 
the line of sight. The direction of the velocity gradient for the velocity maps in Fig. [9] and 
Fig. [To] is shown by the line enclosed in the circle in the lowest left corner of each map. The 
projected specific angular momentum is given by j2D = Jc,proj/Mc, where Jc,proj and are 
the estimated projected angular momentum, and mass of the core, respectively. By making 
the assumption of uniform rotation, as is commonly done in the observations, one can write 
that Jc,2D = I ^, where I = p Mc R^2D is the moment of inertia of the core, Rc,2D the value 
of the projected radius of the core taken to be the square root of its projected surface, and 
p is a factor of order unity which depends on the morphology and radial density profile of 
the core. The angular velocity of the core is then estimated as being f2 = Q/smi, where i is 
the inclination of the rotation axis to the line of sight. Thus j2D is given by: 



32D = P Rc,2D —■ (9) 
sm 1 

An assumption about the inclination of the cores and the value of sin i has to be made. 
We follow Caselli et al. 2002 and use sin i = tt/A, which is the average value of sin i over 
all possible inclinations (Goodman et al. 1993 used sini = 1. The result of this is just a 
systematic decrease of the estimated j2D by a factor of 4/7r ^ 1.27). In the observations, the 
assumption is also often made that cores have a uniform density which implies that p = 2/5. 
Using these two approximations, j2D becomes: 



J2D = ^Rl2D Q- (10) 

5 vr 

Fig Hi] displays the ratio of the 3D and 2D estimates of the specific angular momentum 
Cj = izo/ 32D for projections of the C160 cores along the three axis of the box at t = 1.152 
Myrs = 0.417 tff^ci and t = 1.170 Myrs = 0.424 tff^ci for simulation Bl (69 cores and 207 
velocity maps in total for the two combined time-steps, left figure) and at t = 1.159 Myrs 
= 0.420 t//,c/ for simulation B2 (46 cores and 138 velocity maps in total, right figure). In order 
to avoid being affected by numerical noise, we have used C160 cores that have a minimum of 
5 grid cells (on the 1024^ grid) in any given direction. From Fig. [TTl it is possible to observe 
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that the value of Cj is systematically < 1 for all cores and projections. The Cj distribution 
peaks at ~ 0.1 in simulation B2 and displays an extended peak in the region 0.03 — 0.2 in 
simulation Bl. A two sided Kolmogorov-Smirnov test indicates that there is a probability of 
38 percent for the two data sets in fig. [TT]of being drawn from the same distribution, and 87 
percent when the data is restricted to the range Cj < 1. In terms of the median values, the 
median value of Cj for the C160 cores in simulation Bl for the epoch displayed in Fig. [TTlis 
0.28 ± 0.13 and the median value of the C160 cores in simulation B2 is 0.22 ± 0.11. 

Whether the correction factor Cj from j2D to j^o is considered in terms of its charac- 
teristic value (~ 0.1) or in terms of its median value (~ 0.25), the aim of Fig. [11] is to show 
that the observed shift between the observed and numerically derived distributions of the 
specific angular momentum as seen in Fig. O is most likely caused by the calculation method 
of the specific angular momentum using the global gradient method under the assumption 
of uniform rotation of the cores. We thus conclude, based on Fig. [5] and Fig. [11] that the 
observational determinations of the specific angular momentum tend to overestimate the 
true value of the specific angular momentum by a factor of at least 4 — 5 but most likely 
produces an overestimate of j by a factor of 8 — 10. The origin of the 2D-3D discrepancy 
stems from the fact that the observational method is based on the global gradient method 
assuming uniform rotation, whereas the measurement of the angular momentum in the in- 
trinsic space is a summation running over all parcels of gas in the cores with their more 
complex dynamical behavior. 



4.2. The size-specific angular momentum relation 

A simplified method in theoretical works of assessing the specific angular momentum 
of cores is through the use of a size-specific angular momentum relation [j^o oc R^). If 
motions in the cores that are assimilated to rotational motions are purely due to turbulent 
motions following a Larson-like velocity dispersion-size relation of the form cTc oc R^, then, the 
angular velocity Q = <Jc/Rc will be given hjQcc R^~^ and the specific angular momentum by 
jsD oc R].'^^ . For (3 = 0.38 as observed initially by Larson (1981), the expected dependence 
of is oc Rl'^^. On the other extreme hand, a rigid-body rotation (i.e., Q is constant) 
implies that j^n oc R^. 

Such relations have been derived observationally by several groups (e.g.. Goldsmith & 
Arquilla 1985; Goodman et al. 1993; Phillips 1999). Bodenheimer (1995) combined the data 
of Goldsmith & Arquilla (1985) and Goodman et al. (1993) and found that the exponent 
of the size-specific angular momentum relation is of the order of A ~ 1.6. Phillips (1999) 
compiled a large number of published molecular cloud data and their substructure of clumps 
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and cores and found A to be of the order ~ 1.43 and close to 0.96 for flattened systems. 
Fig. [12] and Fig. [13] display, at a few selected epochs, the intrinsic specific angular momenta 
j3£) of the cores as a function of their characteristic size Rc for C20 cores (top) and C160 
cores (bottom) in simulations Bl and B2, respectively. We perform a robust least absolute 
deviation fit to the data displayed in Fig. [12] and Fig. [13] in addition to the same data at 
the additional epochs we analyze in order to derive the exponent and normalization of the 
relation j^D = A R'^. Fig. [H] displays the time evolution of A and A for the two populations 
of cores in the two simulations. Both A and A are observed to decrease as a function of 
time and most importantly so for the C160 cores. The value of A is observed to decrease 
from around ~ 2.5 down to ~ 1.8 on a timescale of the order of ~ O.btff^ci- This decrease 
in the specific angular momentum indicates a loss of the momentum as the cores evolve and 
become, on average, more gravitationally bound. 

Following our earlier approach of generating observational counterparts to the relations 
studied in the intrinsic three-dimensional space, we construct the same R2D — J2D relation 
of the cores, where R2D is the size of the core seen in projection and taken to be the square 
root of its surface and j2D is the same, earlier defined, specific angular momentum derived 
from the velocity maps of the cores. Fig. [TB] displays the R2D — 32D relation for the C160 
cores at nearly similar epochs in the two simulations Bl(left) and B2 (right). In addition 
to the higher normalization discussed in the previous section, the slopes of the 320 — Rc,2D 
relation are found to be shallower than the corresponding j^£, — R^. relations observed for the 
same cores at the same epochs. Whereas the exponent A of the js^ — Rc is found to be equal 
to A ~ 1.8 and A ~ 2 at these epochs in simulations Bl and B2 respectively, the exponent 
of the i2D — Rc,2D relation, derived from a robust least absolute deviation fit to the data in 
Fig [151 are Asd = 1-06 ± 0.28 and X2D = 1-32 ± 0.21 (these values are X2D = 1-39 ± 0.49 and 
A2D = — 1.28±0.11, respectively when derived a minimization of the Chi-square error statistic 
method). This shows that the observational bias in the overestimate of the specific angular 
moment discussed in § 14. II affects in a more significant way the smallest cores and this results 
in a value of X2D smaller than X^o- Finally, these results suggest that if a simple estimate 
of the specific angular momentum is to be made using a size-specific angular momentum 
relation, a value of A ~ 2 should be used for the exponent of this relation rather than the 
smaller values derived in the observations. The derived values of A in the three-dimensional 
space which yields values of A ~ 1.8 — 2 should nevertheless not be interpreted as the cores 
being in a state of rigid body-rotation as this is not supported by the observations of their 
projected velocity maps. 
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5. CONCLUSIONS 

In this paper, we have analyzed the rotational properties of dense molecular cloud 
cores formed in two magnetized, self-gravitating molecular cloud simulations with a decaying 
turbulence. The two simulations differ by the strength of the magnetic field in the clouds with 
one cloud being mildly magnetically supercritical and the other being strongly magnetically 
supercritical. Our results show that the formation efficiency of dense cores is strongly reduced 
with increasing importance of the magnetic field in the cloud (going down from 33 percent 
per free-fall time in the strongly supercritical cloud to 6 percent for the mildly supercritical 
cloud). We also observe that the median value of the specific angular momentum of the 
high density cores in the mildly supercritical simulation is smaller than the values derived 
for cores in the strongly supercritical simulation. This result is consistent with the fact that 
magnetic braking which leads to angular momentum loss is playing a more important role 
in the cloud where the magnetic field is stronger. 

We have focused our attention on the discrepancies that may arise between estimates 
of the specific angular momentum of the cores derived from the global velocity gradient 
method commonly used in the observations and its the true value measured in the intrinsic 
three-dimensional space. In order to derive the specific angular momentum of the cores fol- 
lowing the observational procedure, we generate synthetic velocity maps of the cores along 
three different projections. The global velocity gradient of the cores is measured from the 
velocity map using the VFIT routine employed initially by Goodman et al. (1993). The 
specific angular momentum is then calculated using the global velocity gradient value under 
the assumption of uniform rotation of the cores. We find, in the two simulations, that the 
distributions of the ratio of the specific angular momentum determined in the intrinsic 3D 
space to the one derived from projected velocity maps peaks at values around ~ 0.1. This 
may well explain the difference by a factor ~ 10 that is observed between the distribution of 
specific angular momentum derived from the intrinsic data in our simulations and the corre- 
sponding real observations using the NH3 and N2H+ molecules of roughly similar excitation 
density than the density thresholds used to identify the cores in the simulations. 

We suggest that the origin of this discrepency (between 2D and 3D) hes in the fact 
that contrary to the intrinsic determination of j which sums up the individual gas parcels 
contributions to the angular momentum, the observational determination of j is based on a 
measurement on the global velocity gradient under the hypothesis of uniform rotation which 
smoothes out the complex fluctuations present in the three-dimensional velocity field. We 
therefore suggest that previous measurements of the specific angular momentum of the cores 
overestimate its true value and that a correction factor of ~ 10 should be applied to these 
measurements as well as to new determinations of the specific angular momentum when 
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using the global gradient method adopted so far in the observations. As already stressed 
by other groups (e.g., Padoan et al. 1998,2000; Pichardo ct al. 2000; Ostriker et al. 2001; 
Ballesteros-Paredes & Mac Low 2002; Pineda et al. 2009; Federrath et al. 2010; Shetty et 
al. 2010) our work further highlights the importance of generating synthetic observations 
from three-dimensional numerical simulations that can be compared to real observations. 
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Table 1. Properties of the populations of C20 cores at two selected epochs in simulation 
Bl. From left to right, the columns represent: the timesteps (1 t//,cZ = 2.76 Myrs), the 

sonic Mach number A4, the number of cores Nc, the maximum number density in the box 
in units of the average number density, the median mass, the median size, the median 
specific angular momentum, and the median rotational parameter. 



tiniestep 


M 


Nc 


leaver 


Mc,„ed (Mo) 


Rc.med (pc) 




Logif){Prot,med) 


t = 0.293 tff^ci 
t = 0.424 tff^^i 


2.97 
3.36 


634 
530 


1904 
359319 


0.0060 ±0.0041 
0.0045 ± 0.0034 


0.021 ±0.011 

0.019 ±0.010 


(2.52 ± 1.18) X 10^" 
(1.76 ±0.81) X 10^0 


-1.59 ±0.09 
-1.69 ±0.10 



Table 2. Similar to Tab. [H but for the C20 cores of simulation B2. 



tiniestep 


M 




Umax 


Mc^med (Mq) 


RcATied (pc) 


jao^med (cni^ s 1) 


Logio{f3rot,med) 


t = 0.280 tff^ci 
t = 0.420 tff^ci 


3.37 
2.59 


640 
451 


2.68 X 10^ 
6.26 X 10^ 


0.0053 ±0.0041 
0.0070 ± 0.0052 


0.020 ±0.010 
0.022 ±0.011 


(2.02 ±0.99) X 10^0 
(2.61 ±1.09) X 10^0 


-1.62 ±0.14 
-1.85 ±0.08 
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Table 3. Similar to Tab. [H but for the C160 core of simulation Bl. 



timestep 


M 


Nc 


naver 


Me,™ed (Mq) 


Rc.med (pc) 


hD.med [cm? S ^) 


Logi(){(3rot,med) 


t = 0.293 tff^ci 
t = 0.424 tff^ci 


2.97 
3.36 


166 
170 


1.90 X 10^ 
3.59 X 10^ 


0.0061 ± 0.0050 
0.0096 ± 0.0009 


0.010 ±0.004 

0.012 ±0.006 


(8.76 ±4.11) X 10^9 
(14.4 ±6.33) X 10^9 


-2.17±0.13 
-2.19 ±0.14 



Table 4. Similar to Tab. [H but for the C160 cores of simulation B2. 



timestep 


M 


Nc 


nmax 
leaver 


Mc^med (Mq) 


Rc,med (pc) 


jsD.med (cm^ S ^) 


Logio{f3rot,med) 


t = 0.280 tff,,i 
t = 0.420 tff^ci 


3.37 
2.59 


276 
171 


2.68 X 10^ 
6.26 X 10^ 


0.009 ±0.007 
0.015 ±0.001 


0.011 ±0.005 

0.014 ±0.007 


(12.3 ±5.53) X IQis 
(14.5 ±6.11) X 10^^ 


-2.23 ±0.09 
-2.46 ±0.07 
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Fig. 1. — Time evolution of the number of cores identified at different density tliresholds in 
the molcular cloud models Bl(left, model with (3p = 0.1) and B2 (right, model with (3p = 1). 
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Fig. 2. — Time evolution of the fraction of the simulation box mass that is found into dense 
cores identified at the thresholds of Uthr = 20 Uaver (C20 cores) and 160 Uaver (C160 cores) 
in the two simulations Bl and B2. 
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Fig. 5. — Distributions of the specific angular momentum of the C20 and C160 cores at five 
different epochs in simulations Bl (left) and B2 (right). The NH3 cores observations are 
those obtained by Goodman et al. (1993) and Barranco & Goodman (1998), and the N2H"'" 
ones are those obtained by Caselli et al. (2002a). 
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Fig. 6. — Distributions of the f3rot values of the C20 and C160 cores at five different epochs 
in simulation Bl (left) and B2 (right). The NH3 cores observations are those obtained by 
Goodman et al. (1993) and Barranco & Goodman (1998), and the N2H"'" ones are those 
obtained by Caselli et al. (2002a). 



0.25 



0.30 



- 3.0 - 

in units of tf 
0.35 



0.40 



0.45 



0.50 



- ^P=0.1 

10^1 - nthr=20 n,,,, 



NH3 cores 



^20 



H \ H 








10 



21 



NjH cores 



10 



20 




)(- — 



_ — ^( 



0.6 



1.0 

time (in Myrs) 



1 .2 



1 .4 



Fig. 7. — Time evolution of the median value of the specific angular momentum for the C20 
cores (top) and C160 cores (bottom) in simulations Bl and B2. In the top panel the median 
value of the specific angular momentum of the NH3 of Goodman et al. (1993) and Barranco 
& Goodman (1998) is over-plotted to the data. In the bottom panel the median value of the 
specific angular momentum of the N2H"'" cores of CaseUi et al. (2002a) is over-plotted to the 
data. 




Fig. 8. — Time evolution of the median value of the rotational parameter for the C20 (top) 
and C160 cores (bottom) in simulations Bl and B2. In the top panel the median value of the 
rotational parameter of the NH3 cores of Goodman et al. (1993) and Barranco & Goodman 
(1998) is over-plotted to the data. In the bottom panel the median value of the rotational 
parameter of the N2H"'" cores of CaseUi et al. (2002a) is over-plotted to the data. 
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Fig. 11. — Distribution of the ratio of the intrinsic specific angular momentum J^d to the 
specific angular momentum derived following the method adopted in the observations j2D 
for C160 cores identified at nearly similar epochs in simulations Bl (left) and B2 (right). 
The Bl data combines cores identified at t = 1.152 Myrs = 0.417 t//,d and t = 1.170 
Myrs = 0.424 tff^ci (in total, 69 cores and 207 velocity maps) and the B2 data displays the 
distribution for cores identified at t = 1.159 Myrs = 0.420 tff^ci (in total, 46 cores and 138 
velocity maps). 
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Fig. 14. — Time evolution of the exponent and logarithm of the coefficient of the Log^j^i)) = 
A + X Log{Rc) relation for C20 and C160 cores identified in simulations Bl and B2. Time 
is given in physical units in the lower horizontal axis and in units of the free-fall time of the 
cloud in the upper horizontal axis. 
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Fig. 15. — Projected radius i?c,2D-projected specific angular momentum j2D relation for C160 
cores identified at nearly similar epochs in simulations Bl (left) and B2 (right). 



